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Abstract 

The quark-antiquark gauge invariant Green function is studied through its dependence 
on Wilson loops. The latter are saturated, in the laige-Nc limit and for large contours, 
by minimal surfaces. A covariant bound state equation is derived which in the center- 
of-mass frame and at equal-times takes the form of a Breit-Salpeter type equation. The 
large- distance interaction potentials reduce in the static case to a confining linear vector 
potential. In general, the interaction potentials involve contributions having the structure 
of flux tube like terms. 
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The Wilson loop PP is defined as the trace in color space of the path-ordered phase 
factor of the gluon field on a closed contour C : 



1 -iq <p dz'^AJx) 

$(C) = — trePe Jc '^^ \ (1) 

Its vacuum expectation value, denoted W{C), 

W{C) = mC)), (2) 

is a functional of the contour C. Loop equations were obtained and studied by Polyakov 
j2] and Makeenko and Migdal [3 IH E]- The Wilson loop essentially satisfies two types 
of equation, which are equivalent to the QCD equations of motion: The Bianchi identity 
and the loop equations (or Makeenko-Migdal equations). A third property, factorization, 
is obtained in the \arge-Nc limit P for two disjoint contours: W{Ci, C2) = W{Ci)W{C2)- 
Further simplification is obtained in the large- A^^c limit of the theory. In that limit, 
for large contours, i.e., at large distances, nonperturbative asymptotic solutions to the 
Wilson loops are represented by the minimal surfaces having as supports the loop contours 
[21 [7j. Therefore, if one is interested only in the large-distance properties of the theory, 
saturation of the Wilson loop averages by minimal surfaces provides a correct description 
of the theory in this regime. In that case, the Wilson loop average can be represented by 
the following functional of the contour C : 

=e-^^^(^), (3) 

where a is the string tension and A{C) the minimal area with contour C. 

Minimal surfaces also appear as natural solutions to the Wilson loop averages in two- 
dimensional gauge theories [S]. 

To deal with the quarkonium bound state problem, one starts with the two-particle 
gauge invariant Green function for quarks qi and q2 with different flavors and with masses 
mi and 7712: 

G{XI,X2]X[,X'2) = {^p2{x2)U{x2,Xi)i)i{xi)i)^{x\)U{x[,x'2)'llj2{,x'2))A,qi,q2- (4) 

Here, t/(x2,xi) is the path-ordered phase factor, 

-ig j ^ dz^A^{z) 
U{x2, xi) = Pe "'2:1 ^ (^5) 

taken along the straight-line X1X2 (and similarly for U {x'i,X2)). Integrating in the large- iVc 
limit with respect to the quark fields, one obtains: 

G{Xi,X2;x[,X2) = -{tlcU{x2,Xi)Si{A;Xi,x[)U{x[,X2)S2{A;X2,X2))A, (6) 
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where Si{A) and 5*2 (A) are the quark and antiquark propagators in the presence of the 
external gluon field and trc designates the trace with respect to the color group. The 
quark propagator S{A) satisfies the equation 

(i'y.dx — m — g'y.A{x)^S{A; x, x') = i6^{x — x). (7) 

At this stage, the Green function G can schematically be represented as in Fig. ^ 

xi Si (A) x'l 




Figure 1: Schematic representation of the two-particle Green function. 



In order to make the Wilson loop structure of G apparent, we adopt for the quark 
propagator in the external gluon field a representation based on an explicit use of the 
phase factor along straight lines [7]. Introducing the gauge covariant composite object 
S{A]x,x'), made of a free fermion propagator So{x — x') (without color group content) 
multiplied by the path-ordered phase factor U (x, x') [Eq. (0)] taken along the straight 
segment x'x, 

S{A-x,x') = Sq{x-x')U{x,x'), (8) 

one shows that the quark propagator S{A\x,x') in the external gluon field satisfies the 
following functional integral equation in terms of S: 

S(A- X, x') = Si A; X, x') - I d^x"S{A- x, x")7" C dX (1 - A)— 4tt5(A; x", x'), (9) 

J Jo ox°'{Xj 

where the segment x"x' has been parametrized with the parameter A as x(A) = (1 — A)x" + 
Xx' and where the operator 6/6x°'{X) acts on the factor U of S, along the internal part of 
the segment x"x', with x' held fixed. That equation is diagrammatically represented in 
Fig.El 

A similar equation in which the roles of x and x' are interchanged also holds. Those 
equations lead to iteration series for S in which the gauge covariance property is main- 
tained at each order of the iteration. 

Use of the above representations for the quark propagators in Eq. (jnj leads for the 
two-particle Green function to a series expansion where each term contains a Wilson loop 
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S(A) S{A) 
= + 

Figure 2: Diagrammatic representation of the integral equation satisfied by the quark 
propagator in the external gluon field. The cross represents the action of the functional 
derivative 6/6x{\). 




along a skew-polygon: 

oo 

G = E Gij, (10) 

where Gij represents the contribution of the term of the series having {i — 1) points of 
integration between xi and x'l {i segments) and (j — 1) points of integration between X2 
and X2 (j segments). We designate by Cij the contour associated with the term Gij. A 
typical configuration for the contour of 6*4 3 is represented in Fig. |21 




Each segment of the quark lines supports a free quark propagator and except for the 
first segments (or the last ones, depending on the representation that is used) the Wilson 
loop is submitted to one functional derivative on each such segment. One then uses for 
the averages of the Wilson loops appearing in the above series the representation with 
minimal surfaces [Eq. Q]. 

The Green function G satisfies the following equation with respect to the Dirac oper- 
ator of particle 1 acting on xi: 

(«7-%i) - mi)G{xi,X2;x[,X2) = -i{ticU{x2, Xi)S^{xi - x[)U{x[,X2)S2{x'2,X2))a 
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-z7"(tre/^'cia(l-a)^^^^5i(a;i,a;;)t/(x;,4)^2(4,^2))A, (11) 

where the segment X1X2 has been parametrized with the parameter a as x{a) = (1 — 
a)xi + ax2', furthermore, the operator d/6x°' does not act on the exphcit boundary point 
Xi of the segment, this contribution having been cancelled by the contribution of the 
gluon field A coming from the quark propagator Si. A similar equation also holds with 
the Dirac operator of particle 2 acting on X2- Representation Q for the quark propagator 
can then be used in the above equation and its partner satisfied by the two-particle Green 
function G. One obtains two compatible equations for G where the right-hand sides 
involve the series of the terms Gij of Eq. (jlUj) and their functional derivative along the 
segment X1X2. In order to obtain bound state equations, it is necessary to reconstruct 
in the right-hand sides the bound state poles contained in G p. In a;-space, bound 
states are reached by taking the large separation time limit between the pair of points 
{xi,X2) and {x[,X2) [IDl- To produce a bound state pole, it is necessary that there be a 
coherent sum of contributions coming from each Gij, since the latter, taken individually, 
do not have poles. Each Gij involves a corresponding minimal surface Aij on the contour 
of which act various functional derivatives. Those can be classified according to their 
possible irreducibility properties. Reducible contributions are those which are parts of 
the definition of the series of G. It does not seem possible to sum all these terms to 
reproduce exactly G with some kernel acting on it in the right-hand sides. However, for 
large separation time limits one can isolate terms that contribute to the pole terms. One 
notices that the derivative along the segment X1X2 acts on areas Aij with contour Gij 
which are different from one term of the series to the other (the number of segments being 
different). To have a coherent sum of those contributions it is necessary to expand each 
such derivative term around the derivative of the lowest-order contour Ci^i, represented 
in Fig. m 




X2 

Figure 4: The lowest-order contour Ci 1 and its minimal surface 1. 
It is that term that can be factorized and can lead through the summation of the 
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factored series to the reappearance of the Green function G and to its poles. The re- 
maining terms do not lead to pole terms. Similarly, two derivative contributions should 
be expanded around the lowest-order contribution coming from the contours C2,i or Ci^2, 
and so forth. 

In general, the derivative of the areas along X1X2 depends among others on the slope 
of the areas in the orthogonal direction to xiX2- One then associates that slope with 
the quark momenta. Taking then the large separation time limit and equal times in the 
center-of-mass frame, one ends up with a covariant three-dimensional equation, having 
the structure of a Breit-Salpeter type equation [TTJ [T2] and where the interaction kernels 
or potentials are given by various functional derivatives involving at least one derivative 
along the segment XiX2- Keeping for the potentials the terms containing one functional 
derivative of the area Ai i [Fig. the equation takes the form [7j 

Po - {hio + h2o) - 7io7i - 72072 ^2/.] ^(x) = 0, (12) 

where is a 4 x 4 matrix wave function of the relative coordinate x = X2—X1 considered at 
equal times, Pq the center-of-mass total energy and hio and /z2o the quark and antiquark 
Dirac hamiltonians; the Dirac matrices of the quark (with index 1) act on from the left, 
while the Dirac matrices of the antiquark (with index 2) act on ip from the right. The 
potentials Ai and A2 are defined through the equations 

-.')^. 4.,=./;d.v'^. (13) 

x(cr') belonging to the segment XiX2- 

The time components of Ai and A2 add up in the wave equation. For their sum, one 
has the expression (in the cm. frame) 

El Eo ( / E-\ , , E2 




+ (e(Pio) - e(P2o))' ^^^^ 

(14) 

Here, r = \f^^ E^ = ym^ + p^^ a = 1,2, with the quark masses, p the cm. momen- 
tum, p = (p2 — Pi)/2, L the cm. orbital angular momentum, and e(pio) and e(p2o) the 
energy sign operators of the free quark and the antiquark, respectively: 

e{Pao) = ^, a =1,2. (15) 
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The space components of Ai and A2 are orthogonal to x. The expression of Ai is (in 
the cm. frame): 



Ai = —err 




+ ^2 V r^El 
Here, p* is the transverse part of p with respect to x: 



(16) 



P = p-x— x.p. (17) 
x^ 

The expression of A2 is obtained from that of Ai by an interchange in the latter of the 
indices 1 and 2 and a change of sign of p*. 

For sectors of quantum numbers where = 0, the expressions of the potentials 
become: 

^10 + A20 = ^(e(pio) + e(P2o))crr, (18) 
1 /I „ . 1 



M = +^^{l{Ei + E,)-'^E2)p'ar. (19) 

The potentials are generally momentum dependent operators and necessitate an ap- 
propriate ordering of terms. 

From the structure of the wave equation (fT^ and the expressions of the potentials, one 
deduces that the interaction is confining and of the vector type. However, compared to the 
conventional timelike vector potential, it has additional pieces of terms contributing to the 
orbital angular momentum dependent parts. A closer analysis of those terms shows that 
they can be interpreted as being originated from the moments of inertia of the segment 
X1X2 carrying a constant linear energy density equal to the string tension. The interaction 
potentials are therefore provided by the energy-momentum vector of the segment joining 
the quark to the antiquark, in similarity with the color flux tube picture of confinement. 
An analogous equation had also been proposed by Olsson et al. on the basis of a model 
where the quarks are attached at the ends of a straight string or a color flux tube ITl] . 
A similar conclusion had also been reached by Brambilla, Prosperi et al. on the basis of 
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the analysis of the relativistic corrections to the nonrelativistic hmit of the Wilson loop 

For heavy quarks, one can expand equation (fT^ around the nonrelativistic limit and 
obtain the hamiltonian to order [Zj. 

The relativistic corrections to the interquark potential arising from the Wilson loop 
were analyzed and evaluated in the literature by Eichten and Feinberg ^H], Gromes [T^ . 
Brambilla, Prosperi et al. [T^ 1161 IT7j . Brambilla, Pineda, Soto and Vairo PU] . 

The Wilson loop approach was also used for the study of quarkonium systems by 
Dosch, Simonov et al. with the use of the stochastic vacuum model 

In conclusion, the saturation of the Wilson loop averages in the large-A^c limit by 
minimal surfaces provides a systematic tool for investigating the large- distance dynamics 
of quark-antiquark bound state systems. 
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